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Abstract 

One of the main drivers for new waveforms in future 5G wireless communication systems is to handle efficiently the 
variety of traffic types and requirements. In this paper, we introduce a new random access within the standard acquisition 
procedures to support sporadic traffic as an enabler of the Internet of Things (loT). The major challenge hereby is to cope 
with the highly asynchronous access of different devices and to allow transmission of control signaling and payload ”in 
one shot”. We address this challenge by using a waveform design approach based on bi-orthogonal frequency division 
multiplexing where transmit orthogonality is replaced in favor of better temporal and spectral properties. We show that 
this approach allows data transmission in frequencies that otherwise have to remain unused. More precisely, we utilize 
frequencies previously used as guard bands, located towards the standard synchronous communication pipes as well as in 
between the typically small amount of resources used by each loT device. We demonstrate the superiority of this waveform 
approach over the conventional random access using a novel mathematical approach and numerical experiments. 
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I. Introduction 


The Internet of Things (loT) is expected to foster the development of 5G wireless networks and requires efficient 
access of sporadic traffic generating devices. Such devices are most of the time inactive but regularly access the Internet 
for minor/incremental updates with no human interaction, e.g., machine-type-communication (MTC). Sporadic traffic 
will dramatically increase in the 5G market and, obviously, such traffic should not be forced to be integrated into the 
bulky synchronization procedure of current 4G cellular systems m, 0. 

The new conceptional approach in this paper is to use an extended physical layer random access channel (PRACH), 
which achieves device acquisition and (possibly small) payload transmission ”in one shot”. Similar to the implementation 
in UMTS, the goal is to transmit small user data packets using the PRACH, without maintaining a continuous connection. 
So far, this is not possible in LTE, where data is only carried using the physical uplink shared channel (PUSCH) so that 
the resulting control signaling effort renders scalable sporadic traffic (e.g., several hundred nodes in the cell) infeasible. 
By contrast, in our design a data section is introduced between synchronous PUSCH and standard PRACH, called D- 
PRACH (Data PRACH) supporting asynchronous data transmission HI . Clearly, by doing so, sporadic traffic is removed 
from standard uplink data pipes resulting in drastically reduced signaling overhead and complexity. In addition, this 
would improve operational capabilities and network performance as well as user experience and life time of autonomous 
MTC nodes 0, 0. Waveform design in this context is a very timely and important topic 0, 0, o. Of particular 
importance is also the line of work in the EU projects METIS (www.metis2020.eu) and 5GNOW (www.5gnow.eu). 

We assume that each D-PRACH’s data resource contains only a very few number of subcaniers (about 5-20 subcarri¬ 
ers). In addition, in a 5G system, we can expect that there is a massive number of MTC devices, which will concurrently 
employ these data resources in an uncoordinated fashion. In the simplest approach, the D-PRACH uses the guard bands 
between PRACH and PUSCH, which is the focus of this papeij^ We show that waveform design in such a setting is 
necessary since the OEDM waveform used in LTE cannot handle the highly asynchronous access of different devices 
with possible negative delays or delays beyond the cyclic prefix (CP). Clearly, guards could be introduced between 
the individual (small) data sections and to the PUSCH which, though, makes the approach again very inefficient. Our 
results indeed show that up to four subcarriers can be obtained compared to a standard 4G OFDM setting. 

Eor waveform design, we propose a bi-orthogonal frequency division multiplexing (BEDM) based approach where 
we replace orthogonality of the set of transmit and receive pulses with bi-orthogonality, i.e., they are pairwise (not 
individually) orthogonal. Thus, there is more flexibility in designing the transmit prototype pulse (or waveform), e.g., in 
terms of side-lobe suppression and robustness to time and frequency asynchronisms 0. The BEDM approach together 
with a suitable waveform is well suited to sporadic traffic, since the PRACH symbols are relatively long so that 
transmission is very robust to (even negative) time offsets. In addition, BEDM is also more robust to frequency offsets 
in the transmission, which, as well-known, typically sets a limit to the symbol duration in OEDM transmission. Einally, 
the concatenation of BEDM and OEDM symbols together in a frame requires a good tail behavior of the transmit pulse 
in order to keep the distortion to the payload carrying subcarriers in PUSCH small. The excellent and controllable 
tradeoff between performance degradation due to time and frequency offsets is the main advantage of BEDM with 

^Notably, in an extended setting this region can be enlarged (by higher layer parameters) but, clearly, at some point new efficient channel estimation 
must be devised different to the proposal in this paper. Recent results in METIS and 5GNOW have outlined a sparse signal processing approach to 
cope with this situation 0, E). 
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Fig. 1. PRACH (blue) and PUSCH (red) regions. A guard interval (GI) separates PUSCH from PRACH in LTE (gray). Parts of this area are used 
for data transmissions of asynchronous users (indicated in green) in a novel D-PRACH, whose size can be variably determined by MAC. Rows in 
this illustration do not represent tme subcamer quantities. 


respect to conventional OFDM. 

We investigate the performance of the proposed approach using mathematical analysis and numerical experiments 
where, for comparison, a standard LTE system serves as a baseline. A part of the numerical results were already 
presented in m . We show how the new approach can actually reduce the interference to the PUSCH region, experienced 
in particular by users whose resources are close to the PRACH. Moreover, we demonstrate that the performance in the 
new D-PRACH region is signihcantly improved by the pulse-shaping approach when multiple, completely asynchronous, 
users transmit data in adjacent frequency bands. 

A. System model 

We consider a simple uplink model of a single cell network, where each mobile station and the base station are 
equipped with a single antenna. We assume there exist two channels -in LTE terminology- the PUSCH and the PRACH. 
In the PUSCH, the data bearing signals are transmitted from synchronized users to the base station using SC-EDMA. 
A small part of the resources is reserved for PRACH, in which, at the hrst step of the RACH procedure, users send 
preambles that contain unique signatures. In this paper, we mainly deal with the PRACH design, trying to leave PUSCH 
operations as unaffected as possible. Specific system parameters can be found in Table and standard textbooks M- 

The time-frequency resource grid for the described channels is illustrated in Eigure [T] To minimize the interference 
between the channels, several subcarriers on both sides of the PRACH are usually unused and serve as a guard band to 
PUSCH. In this paper, however, the D-PRACH is located here, i.e., some users represented by specihc signatures use 
this region to send data by sharing the small number of available subcarriers. Naturally, these users may be completely 
asynchronous to either PUSCH or to each other, which is a serious challenge for OEDM and will be handled by the 
BEDM approach using the so-called ’’spline” waveform, which has ’’good” localization properties in time and frequency. 

Remark 1. An important property of the system is that each user can use some autonomous timing advance (ATA), 
introduced in ifTOll . with respect to the LTE broadcast signal. This will lead to possibly negative as well as positive time 
delays of each user with respect to the receiver’s ’’reception” window. We will see that this can signihcantly lower the 









































Fig. 2. Illustration of autonomous timing advance (ATA). The spline waveform is in fact the candidate soluti on d iscussed i n this paper and appears 
more robust to both symmetric time and frequency offsets. The curves depicted here are taken from Figure in Section |V-D| and constitute one 
of our main results. 


distortion with the new waveforms since the distortion is ’’shifted” symmetric around the zero, where it is much lower. 
The principle is shown in Figure 

Eventually, it is worth emphasizing that the guard bands in 4G LTE are relatively large so that the application of 
ATA is restricted to relatively demanding settings with large time and frequency shifts. However, future 5G systems 
are expected to have shorter symbol lengths so that the results of this paper are applicable to much less demanding 
scenarios. 


B. Organization 

The paper is organized as follows. In Section we describe the fundamentals of BEDM using Gabor theory and 
introduce the spline pulses. Then, in Section we present our approach to multiuser interference analysis in the 
context of highly asynchronous access and provide examples for the OEDM and spline waveforms. In Section |IV| we 
deal with practical implementation issues for BEDM. In Section |V] we investigate the performance of the proposed 
approach numerically and compare to standard LTE. In Section]^ we summarize our findings and draw some important 


conclusions. 




































11. BFDM System Design 

A. Elements of Gabor signaling 

Conventional OFDM and BFDM can be formulated as a pulse-shaped Gabor multicarrier scheme. For the time- 
frequency multiplexing we will adopt a two-dimensional index notation n = (ni,n 2 ) G Let i denote the imaginary 
unit and p = (/ii,/i 2 ). The baseband transmit signal is then: 

S{t) = '^Xn7nit) = '^Xn{SAn'y){t), ( 1 ) 

nGl n^I 

where (S'^ 7 )(f) := — fii) is a time-frequency shifted version of the transmit pulse 7 , i.e., 7 „ := SauI 

is shifted according to a lattice hi?. The lattice is generated by the 2x2 real generator matrix A and the indices 
n = (ni, 712 ) range over the doubly-countable set I Gl?, referring to the data burst to be transmitted. The coefficients 
Xn are the random complex data symbols at time instant ni and subcarrier index 712 with the property Ejtca;*} = I 
(from now on “ always denotes complex conjugate, •* means conjugate transpose, and x = (... ... )^). We will 

denote the linear time-variant channel by the operator H and by 77 (f) an additive distortion (a realization of a noise 
process). 

In practice, A is usually diagonal, i.e., A = diag(T, F) and the time-frequency sampling density is related to the 
bandwidth efficiency (in complex symbols) of the signaling, i.e., e := | det A“^| = (TF)~^. The received signal is: 

r(t) = {'Hs){t) + n{t) = [ 'E{fi){S^s)(t)d^ + n{t), (2) 

with S : C being a realization of the (causal) channel spreading function of finite support. In the wide-sense 

stationary uncorrelated scattering (WSSUS) assumption ini, the channel statistics is characterized by the second order 
statistics of S, given as the scattering function C : —>■ M+: 

= (3) 

Moreover, we assume E{S(/i)} = 0 and ||C||i = 1. To obtain the (unequalized) symbol Xm on time-frequency slot 
777 G I, the receiver projects on SAm9' 

Xm-={9m,r)= y (A777)i)r(f)(ff, (4) 

using the £2 scalar product (t) By introducing the elements 

dim,n •— {dm: ( 3 ) 

of the channel matrix H G jjje overall transmission can be formulated as a system of linear equations 

X = Hx + n, where n = (..., {gm,n ),... )^ is the vector of the projected noise. We use the AWGN assumption 
such that n is Gaussian random vector with independent components, each having variance := E{|rm, 77 )p}. The 
diagonal elements 

Hm,m = f ( 6 ) 

Here, 

As 7 (m) := { 9 , Sf.'y) = j^g{t) (^^ 7 ) {t)dt (7) 


is the well known cross-ambiguity function of g and 7. 




Example 1. A lattice can be described by a so-called generator matrix, which determines the geometry. For cp-OFDM, 
we can define the matrix 


A = 


T -b Tcp 
0 


0 

T. 


(8) 


At the transmitter, the rectangular pulse 


1 


lit) = Xl-T^[cp],T]it) 


(9) 


is used, with X[-T([cp\,T\ being the characteristic function of the interval [—T[[cp],T]. The receiver obtains the complex 
symbol as 


— 


5(f- 


( 10 ) 


using the rectangular pulse 


git) = ^X[o.T] 


it) 


( 11 ) 


for the removal of the CP. 


B. Completeness, Localization and Gabor Frames 

In the following, we will collect some well-known statements on (bi-infinite) Gabor families Qii, A) := {S'AnTjnGZ^- 
A more detailed discussion of these concepts with respect to multicarrier transmission in doubly-dispersive channels 
can be found for example in ifT^ and lfT3]| . Based on the generator matrix A, we can categorize the following regimes: 
we refer to critical sampling in case det(A) = 1 , while det(A) < 1 and det(A) > 1 is called oversampling and 
undersampling of the time-frequency plane, respectively. Perfect reconstruction for any X C 1?, i.e., {gm,ln) = ^mn 
for all m,n G X?, can be achieved if and only if det(A) > 1. In this case, g is called the dual (bi-orthogonal) pulse 
to 7 with respect to the lattice A. If p = 7 , the Gabor family Qig,A) is an orthogonal basis for its span. However, a 
main consequence from the Balian Low Theorem is, that no well-localized prototype 7 (in both time and frequency) 
can generate a Gabor Riesz basis at det(A) = 1. Thus, any orthogonal or biorthogonal signaling at the critical density 
is ill-conditioned and will be highly sensitive with respect to either time or frequency shifts. Therefore, without further 
constraints on the data symbols, one has to operate in the undersampling regime in a practical scenario. 

Let us now consider the adjoint lattice, generated by A° = det(A)“^ • A, i.e., if ^( 7 , A) refers to undersampling, 
^( 7 , A°) corresponds to oversampling. An important notion here is the concept of a Gabor frame, i.e. ^( 7 ,A°) 
establishes a frame (for L 2 iM.)) if there are frame bounds 0 < < 00 such that: 

All/ll2< (12) 

nGZ2 

for all / G L 2 iR). With 0 < Aj < < 00 we always mean here the best bounds. The Ron-Shen Duality now states 

that the Gabor family C/( 7 , A) is a Riesz basis for its span if and only if ^( 7 , A°) is a frame for L 2 (K)- Furthermore, 
Qig, A°) is a tight frame (A,y = B^f) \f and only if ^( 7 , A) is an orthogonal basis for its span. If < cx) exists in 0 
for the generator matrix A, the sequence of elements in Q ( 7 , A) is called Bessel sequence. A straightforward argument 
shows that B,y is the maximal eigenvalue of the bi-infinite Gram matrix G.y with the components iG,y)m.,n = {irmln)- 
More details on these concepts in multicarrier transmission can be found in iflTll . 







Fig. 3. Illustration of I®* order spline function. It is the result of a convolution of two 0*^ order rectangular' pulses. 


Example 2 ( ifTSl Appendix A]). The Bessel bound of cp-OFDM signaling is given by the operator norm of S.y, which 
is equal to the largest eigenvalue of the Gram matrix i.e., 


where e = 


T^+To 


(C ) -A 

[^y)m,n-0m„u,e f(n2-m2) ’ 

This is a Toeplitz matrix in the frequency slots m2 and 712 generated by: 
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= eM--u;J+l . 


(13) 


( 14 ) 


(15) 


Using this, the Bessel bound for cp-OFDM can be upperbounded by By = 2 . However, of course this is only the 
worst-case estimation. Interestingly, by contrast, the Bessel bound Bg for the receive ’’frame” is exactly unity, i.e.. 


Bg = l. 


C. Spline-based Gabor Signalling 

As mentioned before, the used pulses g and 7 play a key role and should therefore be carefully designed. As we 
consider a BFDM approach, we generate the transmit pulse g according to system requirements and compute the receive 
pulse 7 as the canonical dual (biorthogonal) pulse. For this, we use a method already applied, for example, in M- 
Briefly explained, bi-orthogonality in a stable sense means that g should generate a Gabor Riesz basis and 7 generates 
















the corresponding dual Gabor Riesz basis. From the Ron-Shen duality principle ifTbl follows that 7 has the desired 
property if it generates a Gabor (Weyl-Heisenberg) frame on the so called adjoint time-frequency lattice and that frame 
is dual to the frame generated by g. This can be achieved with the ^'“^-trick explained in El. Side effects such as 
spectral regrowth due to periodic setting when calculating the bi-orthogonal pulses are negligible. 

As a rough and well-known guideline for well-conditioning of this procedure, the ratio of the time and frequency 
pulse widths (variances) at and ct/ should be approximately matched to the time-frequency grid ratio 


T 

F 


(16) 


and this should also be in the order of the channel’s dispersion ratio d. However, since we focus on a design being 
close to the conventional LTE PUSCH and PRACH, of this rule, we consider only ( [T6] l here. 

We propose to construct the pulse g based on the B-splines in the frequency domain, see Figure i3-splines have 
been investigated in the Gabor (Weyl-Heisenberg) setting for example in ifTSl . The main reason for using the B-spline 
pulses is that convolution of such pulses have excellent tail properties with respect to the Li-norm, which is beneficial 
with respect to the overlap of PRACH to the PUSCH symbols. We also believe that they trade off well the time offset 
for the frequency offset performance degradation but this is part of further on-going investigations and beyond the 
conceptional approach here. Because of its fast decay in time, we choose a second order B-spline (the ”tent”-function) 
in frequency domain, given by 


B 2 (/) =5i(/)* Si(/), where (17) 

Si(/):=X[_i.i](/) (18) 

(and * denotes convolution). It has been shown in ifTSl that S 2 (/) generates a Gabor frame for the (a, 6)-grid (translating 
B 2 on aZ and its Fourier transform on bZ) if (due to its compact support) a < 2 and b < 1/2, and fails to be frame in 
the region; 

{a > 2, 6 > 0}U {a > 0,1 < 6 G N}. (19) 


Recall, that by Ron-Shen duality iflhl it follows that the same pulse prototype B 2 {f) generates a Riesz basis on the 
adjoint (|, i)-grid. In our setting, we will effectively translate the frequency domain pulse B 2 {f) by half of its support, 
which corresponds to | = 1, and we will use ^'7 = 1 = 1-25 (see here also Table such that a = |. Therefore, 
it follows that our operation point {a,b) = (|,1) is not in any of two explicit (a, 6)-regions given above. But for 
1.1 < a < 1.9 a further estimate has been computed explicitly for B 2 {f) ifT^ Table 2.3 on p.560], ensuring the Gabor 
frame property up to 6 < 1/a. Finally, we like to mention that for ab < 1/2 the dual prototypes can be expressed again 
as finite linear combinations of B-splines, i.e., explicit formulas exists in this case m. Note that we can choose a 
larger grid in the frequency domain and set & < 1/2 so that a = 8/5, which is a frame (so that the spline is a Riesz 
basis). Hence, due to the spectral efficiency constraint with increasing 1/5, we also decrease the time domain grid such 
that, necessarily, at some point a > 2 (in the dual domain) and so we do not get a Riesz basis (or a frame in the dual 
domain). 

In practice, g has to be of finite duration, i.e., the transmit pulse in time domain will be smoothly truncated: 

/ sin(B 7 rf)' 

V Birt 


git) = 


X[c,d\ {t) ? 


( 20 ) 




Fig. 4. The cross-ambiguity function t) for transmit and receive pulses on the frame length of 4 ms. 


where B is chosen equal to F and parameters c and d align the pulse within the transmission frame. Theoretically, a 
(smooth) truncation in ( |20| would imply again a limitation on the maximal frequency spacing B ll20l . Although the 
finite setting is used in our application, the frame condition (and therefore the Riesz basis condition) is a desired feature 
since it will asymptotically ensure the stability of the computation of the dual pulse 7 and its smoothness properties. 

To observe the pulse’s properties regarding time-frequency distortions, we depict in Figure the discrete cross¬ 
ambiguity function Agbetween pulse g and 7 . 

It can be observed that its value at the neighboring symbol is already far below 10“^. Obviously, the bi-orthogonality 
condition states Ag^j{kT,lF) = SkflStg and ensures perfect symbol recovery in the absence of channel and noise. 
However, the sensibility with respect to time-frequency distortions is related to the slope shape of Ag,.^ around the 
grid points. Depending on the loading strategies for these grid points it is possible to obtain numerically performance 
estimates using, for example, the integration methods presented in lEl. 

Let us introduce a parameter a to scale the width of the spline pulse in frequency. Choosing a large pulse width 
a has the disadvantage of an increasing value of Bg in (41 1 . This is illustrated in Figure]^ It shows the upper frame 
constant Bg for the transmit pulse g vs. its pulse width a. The frame constant is thereby calculated using the LTFAT 
toolbox ED. It can be observed that Bg has its smallest value around a = 0.85, however, it never reaches the lower 


bound Bg = 1. 













Fig. 5. Upper frame constant B^y over a. A minimum of this function can be observed around 1 (more precisely, at a = 0.85), while for increasing 
and, even more, for decreasing values of alpha the value of B-y increases. 


III. Multiuser Interference Analysis 

A. A general approach 

Our system model has to capture that many users each occupy a small number of subcarriers and each of them 
asynchronously (in time and frequency or both) access this resource in an uncoordinated fashion. For a particular time- 
frequency slot m = {mi, m 2 ), we will denote the (random) channel operator as H{m) and the asynchronism as T) (m). 
We assume that H{m) can be estimated using channel estimation procedure while T) (m) cannot. Writing the received 
complex symbol Xm in the absence of additive noise yields 

_ . - . 

~ “f {^m,Tn H m) A ^ ^ ^m.nXn (21) 

n^X ,n^m 

ICI 

where we defined Hm ■= i e., the mean value conditioned on a fixed channel Thus, the transmitted 

symbol Xm will be multiplied by a constant and disturbed by two zero mean random variables (RV), and ICI. The 
first RV Am represents a distortion, which comes from the randomness of T) (m), the second term ICI represents both. 
The mean power of both contributions, conditioned on a fixed channel are := E|^{|AmP} = Pm - 
where Pm ■= E|^{|iTm,mP} and Im ■= E{|ICIp}. Each element of the distortion sum J2n^m ^ni,nXn is given by: 

Hm,n = {gm,P {n)'H{n)yn)- ( 22 ) 


^'As a matter of fact, the expectations depends only on the marginal distribution of 




















Remark 2. Notably, even if T) (n) is the identity (synchronous access) Xm is affected by all other individual contributions 
where the operators depend also on the index n = (^ 1 , 712 ). Hence, the performance for individual slots will be quite 
different, which complicates the situation, and no analytical approach is available so far! If all T) (n) and H{n) are 
independent of n, standard analysis can be used 02 . 

To find a tractable way, we consider the following approach: We assume that the Hm,m can be estimated, and 
we consider the distortion sum Hm,nXn averaged over all the subcarriers m. Obviously, this will average out 

individual interference for a specific subcarrier, but we can assume that these interference terms do not differ much. 
Individual performance is then measured by Hm,m only! Then, we average over the random operators V (n) and 'H{n). 
Let us first consider the sums: 


^ ^ = (23) 

m n^m m n^m 

(24) 

n m m 

< BgEg Y, \\V (n) " E I <3™’ ^ (25) 

n m 

= BgEg Y 112? in) ■ (26) 

n m 


Here, Bg is the Bessel bound of the Gabor family Q{g,A). In the last step, we see that only the ’’action” of the 
operators {I? (n),'H(n)} on 7 is relevant. We have set without loss of generality || 7 ||| = 1 and 1 < llfflli < Eg 
(typically Eg « 1). Next, we compute the expectations and we use = Pm — \Hm\'^' 


^ (27) 

m n^m m 

< BgEg Y E l|2? (n) nin)jjl - ^ E + E ^ 77^ (28) 

n mm 

< BgEg Y E l|2? (rr) - E^(P„ - Dm) (29) 

n m 

= BgEg Y E l|2? (rr) - Y (30) 

n m 

We assume that the asynchronisms cannot increase the received power. For the first term, we estimate 

Y^\\'Ein)n{n)^X < ^E||?7 (n)7n||^ < (31) 

n n n 

according to Q. It remains to bound the second term 

nHm? =X^\n{9m,'D[m)'H{m)^m)\^ (32) 


for some m. For a^b G we abbreviate [a, 5] := 0162 — 02^1 (the symplectic form). Define now the following 

m ) V (to) Sfj,Jm) 

= SXm^ (m) SAmSgj), 


function: 


(33) 

(34) 


which essentially contains the distortion of the /rth contribution in terms of the pulses conjugated by SAm, i e., ’’shifted” 
to TF-slot m in the time-frequency plane. For a fixed channel S, we have Hm = (S, Sm) and on average, with respect 
to H{m), we have; 

E{|if„|2} = (C^,|s^|2). (35) 


Hence, altogether we have the following upper bound on the total expected distortion: 

E ^ (/„ + Dm) < Y, [EgBgWCmWl " {Cm, |s™|")). ( 36 ) 

m£X m£l 

Let us fix the normalization such that Hence, averaging over I, we have proved the following 

theorem. 


Theorem 1. Suppose ||7m||2 = 1 (without loss of generality)^ llffmlli = Eg such that {gm,lm) = 1 (perfect reconstruc¬ 
tion in noiseless case). The average distortion power per subcarrier is upperbounded by 


PI E 'Y. + Dm) f) EgBg 1^1 'Y, {Cm, [Sml ), 


m^X 




(37) 


m^X 


where: 


Sm{-) = ^\H{ 9 m,D (m) 5'(.)7m). 

Example 3. As a special case, assume now a deterministic time-frequency shift T) (m) = S^(m). This distortion is 


non-random and energy preserving, i.e., \\V (m) gm \\2 = IlSmIb- Evaluating the function Sm in (33 i gives: 

Sm{g{m)) = e-^^^^>^^^^’^^Hg,S*AmSuim)SAmS^^m)l) = ’^^Hg , S,^m)S^^m)l) (38) 

^ + S^^m)+^^m)l) (39) 


Hence, in the AWGN case we have: 


I^-I E 3" Dm) f: EgBg Y^ |-^g7(*^(^))l ■ 

' ' mel ' ' mex 


(41) 


B. OFDM 

The cross ambiguity function for 7 and g, as introduced in can be compactly written as (see US) 

~ -^9s(([^l]cpj ^ 2 )), 

where the phase fo = ttvTu is related to our choice of time origin and 


]cp • X [t]cp — \ 0 


r < 0, 

0 < r < r. 


cp. 


T — T T > T 

' ^ cp ‘ ^ Cp' 


(42) 

(43) 


(44) 


The signal quality in the presence of time- and frequency shifts can now be directly obtained from ( |43| ). Apart from 
[•]cp and (the loss in mean signal amplitude due to the CP) (43 1 agrees with the well known auto ambiguity function 


Agg for a rectangular pulse g of width r„. 
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Fig. 6. Simulated interference over frequency offset using different numbers of subcarriers, compared to theoretical bound. In addition, the behavior of 
the two-user case is illustrated, where each user gets half of the available subcarriers and has a different frequency offset. The (aggregate) interference 
converges to the value of the bound when the differences between the two offsets, which are centered around 0.2 • 1250iT2:, and their average become 
smaller. 


Let us also consider time-invariant channels, i.e., (distributional) spreading functions of the form S(p,) = h{iJ,i)S{n 2 ). 
Here, h is called the channel impulse response. If the system exhibits a timing offset i/i with + r^/Jcp = 0 only, the 
time dependency in the cross ambiguity function cancels, thus 




(45) 


and only phase rotations occur (normally corrected by channel estimation and equalization). Contrary to this, time 
offsets with [vi + T^ijcp ^ 0 cause interference. For frequency offsets, interference occurs immediately. For the case 
[i/i -f Tdjcp = 0, the following relation holds: 


'KI>2 


(46) 


Obviously, the frequency offset 1)2 (the time offset i>i) induces a rotating phase over the time slots I (over frequency 
slots k) as one would expect. 

Finally, we need the Bessel bound (by contrast to transmit pulse) for the receive pulse g, which is Bg = 1 (see 
Example 2), and the energy constant, which is Eg = e, so that altogether 

^ W\^(T 4- n \ < F ^ sin^(7ri/™(Tu - |[:/{"]cp|)) 








-X[-T„,T„]([i^r]cp)- 


(47) 
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Fig. 7. Simulated interference over time offsets using different numbers of subcarriers, compared to theoretical bound. 


1) Evaluation: First, we consider only frequency offsets with no additional delay in time. Figure|^shows a comparison 
of this interference bound with simulated curves at different numbers of subcarriers. It can be observed that with an 
increasing number of subcarriers considered, the interference curve gets closer to the theoretical bound. In case of 200, 
or even more 20000, subcarriers, the simulations match the bound almost perfectly. 

The described curves are based on a single frequency offset only, which is the same for all subcarriers. However, Figure 
1^ additionally illustrates the behavior in case of multiple different offsets. For the sake of illustration, we demonstrate 
the case of two offsets here, where each offset applies to an equal share of the available subcarriers. It can be observed 
that with decreasing difference in the offsets, the resulting interference level gets closer to value of the bound at the 
corresponding average of the offsets. 

Let us now consider the case of time delays. Figure compares the interference caused by the asynchronous mode 
of operation to that predicted by the theoretical bound. Again, it can be observed that the numerical results converge 
to the bound when increasing the number of subcarriers. Note that the CP length is 103ms; an smaller offset does not 
produce any interference (however, negative delays do, which is not depicted here). 


C. Spline-based modulation 


Let us now investigate the influence of time and frequency offsets on the spline-based waveform, as carried out for 


OFDM in Section III-B In order to obtain a bound on the influence of time or frequency offsets, we need to evaluate the 


cross-ambiguity function Ag-y. For this, the dual pulse has to be taken into account. However, since this is not desirable. 













Fig. 8. Interference bound for spline waveform based over frequency offset for increasing values of alpha. Note that = 1 is assumed in all 
curves. 


simple estimates of this function are needed. For this purpose, consider the following bound. Define 7 ^ and 

7 ,. := and using ( 5 , 7 ) = 1 and ||p ||2 = Eg-. 

l^g7(M)l = 1 ( 5 ,7m -7 + 7)1 (48) 

= |1 +(5,7m- 7)1 = |1 + (5,7m-7i.+7<.-7)I (49) 

> 1 - 1 ( 5 , 7m-7-^)1 - 1 ( 5 ,7>^ - 7)1 (50) 

> l-£;g||7^-7,||2-F;g||7,-7||2, (51) 

where the RHS constitutes a similarity measure for 7 . Let us take /r = {At, Aui) and 1 / = (0, Aw) and it follows; 


Theorem 2. Suppose ||7m||2 = 1 (without loss of generality), Upmlll = such that {grmlm) = 1 (perfect reconstruc¬ 
tion in noiseless case). Then 

1 ^ 37 ( 5)1 > 1 - II 7 - 7(- - Af )||2 - II 7 - 7(' - Aa;)|| 2 , (52) 

where 7 denotes the Fourier transform of 7. 


From this, we get by standard analysis the following approximation for frequency offsets: 

1/3AwT //\ AujT\ 


Ag^((0,Auj)) > 1- 


27ra 


1 - 


27rQ; / 


(53) 


Using this approximation we can calculate the interference part in ( |4l] i. Figure illustrates the behavior of the spline 
waveform with different pulse widths a. It can be observed, that the interference decreases with increasing pulse width. 
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Fig. 9. Simulated interference and interference bound vs. frequency offsets for the spline waveform. Numerical results for different numbers of 
subcaniers are shown. To foster an easy comparison, the interference bound for OFDM is also depicted. 


Similar to frequency offsets, we can derive: 


Ag^{{At,0))>l- 


2'KaAt 

V^T ■ 


(54) 


1) Evaluation: In Figure]^ we show the simulated interference of the spline waveform together with the corresponding 


bound based on a numerical computation of Ag^ in (41 1 for frequency offsets. In addition, the figure again shows also 
the interference bound for OFDM. While, the results for the spline waveform indicate lower interference than in the 
OFDM case, the bound appears to be less tight, even with large numbers of subcaiTiers. 

Let us now consider the case of time offsets. Figure [T^ depicts the interference bound, again based on a numerical 


calculation of Ag^ in (41 1 and the simulated interference vs. a time offset for the spline waveform. To allow an easy 
comparison, the bound for OFDM is also shown. It should be noted that although the results do not outperform OFDM 
for the positive delays considered in Figure [T0| the behavior is different for negative delays. In Figure [TT] we show the 
influence of negative delays, where the benefits of the spline waveform become obvious. 


IV. Practical Implementation 

In the following, let Tg denote the sampling period, which is equal to l//s, with /g being the sampling frequency. 
In the following discrete model, we let all time indices be multiples of Tg and frequency indices be multiples of F. 
Furthermore, we use N to denote the discrete counterpart of the symbol duration and submit K symbols. Let iVpFT 
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Fig. 10. Simulated interference and interference bound vs. time offset for the spline waveform. Numerical results for different numbers of subcarriers 
are shown. To foster an easy comparison, the interference bound for OFDM is also depicted. 


be the FFT-length in D-PRACH. Note, for some numerical reasons VppT must divide NK. In order to be compliant 
with 4G, we set TF = 1.25. Note that efficient implementations are available in literature ll22ll . 


A. Transmitter 


For the pulse shaped PRACH, additional processing is needed, compared to standard OFDM. In contrast to standard 
processing, we process more than one symbol interval, even if we use only one symbol to carry the preamble. A pulse 
g is used to shape the spectrum of the preamble signal (which is constructed from a Zadoff-Chu (ZC) sequence ||23), 
e.g., to allow the use of PRACH guard bands with acceptable interference. Let P be the length of pulse g. We extend 
the output signal s[n] after the inverse FFT (IFFT) stage by repeating it and taking modulo P to get the same length 
as the pulse g. Given K symbols, each symbol Sfc[n] is pointwise multiplied by the shifted pulse g and superimposed 
by overlap and add, such that we get the baseband pulse shaped D-PRACH transmit signal as: 

K-l 

= X! ^k[n]g[n - kN]. (55) 


In greater detail, this can be written as 


K-l 


= /3 X! X! - kN\t 

k—0 


(56) 


where Xk,i is the Fourier transformed preamble signal (which is constructed from Zadoff-Chu (ZC) root sequences 
in LTE-A) occupying subcaiTiers I G Fi of length Nzc and the actual data occupying subcaiiiers I G F 2 31 the fc-th 













Fig. 11. Simulated interference and interference bound vs. negative and positive time offsets for OFDM and spline waveform. Numerical results 
for different numbers of subcaiTiers are shown. To foster an easy comparison, the interference bound for OFDM is also depicted. 


symbol, and /? is an amplitude scaling factor for normalizing the transmit power. 

B. Receiver 

Assuming an AWGN channel and one user transmitting its preamble signal on the PRACH, the base station obtains 
the superposition of data bearing signals, preamble signal, and noise as 

r[n] = spu W + spR[n] + no[n], (57) 

where spu is the PUSCH data transmit signals, spR is the PRACH preamble transmit signal, and no is Gaussian noise. 

The only difference of our BFDM receiver to the standard OFDM PRACH receiver is the processing before the 
FFT. While in standard PRACH processing the CP has to be removed from the received signal rpR[n], in our BFDM 
receiver, an operation to invert the (transmitter side) pulse shaping is carried out. To be more precise, the K symbols 
of the received signal rpR[n] are pointwise multiplied by the shifted bi-orthogonal pulse 7, such that we have 

[n] = Tk [n]7* [n - kN]. ( 58 ) 

Subsequently, some kind of prealiasing operation is applied to each windowed i.e., 

P / N-pFT — 1 

£ r2[n-lNFFT], 


( 59 ) 









such that we obtain the Fourier transformed preamble sequence at the A:-th symbol and Z-th subcarrier after the FFT 
operation 


TVfft — 1 




"-7 

n ' 


j 27 rnl 


r/ n e . 


(60) 


n—0 


Although we do not employ a CP as in standard PRACH, the time-frequency product of TF = 1.25 allows the signal 
to have guard regions in time and frequency as well. This time-frequency guard regions and the overlapping of the 
pulses evoke the received signal to be cyclostationary which gives the same benefit as the cyclostationarity made 
by CP. Furthermore, it is also shown in ll24ll . that the bi-orthogonality condition of the pulses is sufficient for the 
cyclostationarity and makes it possible to estimate the symbol timing offset from its correlation function. 

C. User Detection 


1) Preamble generation: The preamble is constructed from a ZC sequence as 

. TTumim + 1 ) 


Xu[m\ = exp < -j 


Nzc 


,0 < m < Nzc - 1, 


( 61 ) 


where u is the root index and Nzc is the length of the preamble sequence, which is fixed for all users. Here, we consider 
the case of contention based RACH, where every user wanting to send a preamble chooses a signature randomly from 
the set of available signatures S = {1,..., 64 —A);f}, with N^ being a given number of reserved signatures for contention 
free RACH. Every element of S is assigned to index {u, v), such that the preamble for each user is obtained by cyclic 
shifting the u-th Zadoff-Chu sequence according to Xu^vim] = x„[(m + v Ncs) mod A^zc], where u = 1,..., 
the cyclic shift index and Ncs is the cyclic shift size. Since only V = 


Nzc 

Ncs 


J preambles can be generated from the 


root u, the assignment from S to (u, v) depends on Ncs and on the size of set S. 


2) Signature detection: Given the received signal (57 1 , the PRACH receiver observes the fraction y that lies in the 


PRACH region to obtain the preamble. The receiver stores all available Zadoff-Chu roots as a reference. These root 
sequences are transformed to frequency domain and each of them is multiplied with the received preamble. As discussed 
in Section]^ it approximately holds, as in OFDM, 


Zu[w] = Y[w]X*[w], 


(62) 


where Y[w\ is the received preamble and Xu[w] is the u-th ZC sequence in frequency domain respectively. Using the 
convolution property of the Fourier transform it is easy to show that Zu[w] is equal to the inverse Fourier transform 
of any cross correlation function Zu[d] at lag d. Because the preamble is constructed by cyclic shifting the Zadoff-Chu 
sequence, ideally we can detect the signature by observing a peak from the power delay profile, given by 


Nzc — ^ 

E 

n—O 


■ d]x* [ 


(63) 


Let A^root = be the number of roots that we require to generate 64 — A).f preambles. Then, the signature Si 


and the delay di of user i are obtained hy Si = Vu + 


, (0 < u < A^root) and di = {n mod Ncs) x Ts, 


respectively, where t; is the location of the largest peak in ([6^. 












D. Channel Estimation 


The question remains how to obtain an estimation for the channel also on the new D-PRACH subcarriers. Due to 
our system setup, we assume that the received preamble signal can be written as 

y = D_^-h + e. (64) 

$ 

Thereby, the term e accounts for all interference and noise, D is a diagonal matrix constructed from the coefficients 

of the Fourier transformed preamble and W = F{Ip,Ih) is a sub-matrix of the DFT-matrix F. The set 

Ih := {1, • ■ • ,nh} contains the indices of the first rih columns, and Ip = {ii,... jiN'z.c} contains the indices of the 
central Nzc rows of F. Furthermore, M is the length of the subframe without CP and guard interval, and we assume 
a maximum length Uh of the channel h. 

For simplicity, we consider simple least-squares channel estimation, i.e., we have to solve the estimation (normal) 
equation To handle cases where $ is ill-conditioned, we use Tikhonov regularization. This popular 

method replaces the general problem of minj-ljAa; — y|P by mina;|jAa; — j/|p + |jra:p, with the regularization matrix 
r. In particular, for our model in ( |64| i 

/i = (65) 

is used in place of the pseudo-inverse, where F has to be adapted to the statistical properties of e. We choose F to be 
a multiple of the identity matrix. 

The idea behind the estimation approach is, that the estimated channel is also valid for subcarriers that are adjacent 
to the region that we actually estimate the channel for. Numerical experiments indicate that the MSE is smaller then 
10“^ for up to 200 subcarriers outside of Ip. 

V. Performance Evaluation 

In this section we verify, using numerical experiments, that using the PRACFI guard bands to carry messages is indeed 
practicable. We compare the standard (LTE) PRACH implementation to our proposed spline pulse shaped PRACH. 

A. Simulation Setup 

The simulation parameters, chosen according to LTE specifications, are provided in Table Eor the computation of 7 
we use the LTEAT toolbox, which provides an efficient implementation of the 5'“^-trick 11211 . Due to the properties of 
the pulses, and to fit the strict LTE frequency specification, we allow a small spillover effect from PRACH to PUSCH in 
time. Due to the PRACH pulse length of 4 ms, as depicted in Pigure[2 we simulate the PUSCH over this time interval. 
Eurthermore we use the maximal available LTE bandwidth of 20 MHz. 

In the LTE standard, the power of PRACH is variable and is incrementally increased according to a complicated 
procedure. To allow a meaningful comparison without having to implement to complete PRACH procedure, we choose 
the power of the PRACH such that approximately the same power spectral density as in PUSCH is achieved, as 
depicted in Pigure[^ We simulate multipath channels with a fixed number of three channel taps. Moreover, we assume 
a maximum length of n = 300, which corresponds to a delay spread of roughly 5 /iS, and which implies a maximum 
cell radius of 1.5 km. Eor the transmission in PRACH, we use 4-QAM modulation. Consequently, even in case the 
PRACH power is lower than in Figure [T^ we still have the opportunity to reduce the modulation to BPSK. 


TABLE I 

System SPECincATioN 



PUSCH 

Standard 

PRACH 

pulse shaped 
PRACH 

Bandwidth 

20 MHz 

1.08 MHz 

1.08 MHz 

OFDM symbol duration 

0.67 fis 

SOO fis 

- 

Subcarrier spacing F 

15 kHz 

1.25 kHz 

1.25 kHz 

Sampling frequency fs 

30.72 MHz 

30.72 MHz 

30.72 MHz 

Length of FFT N^ft 

2048 

24576 

24576 

Number of subcarrier L 

1200 

839 

839 

Cyclic prefix length Tqp 

160 Ts 1st 

144 Ts else 

3168 

0 

Guard time Tg 

0 

2976 Ts 

0 

Pulse Length P 

- 

- 

4 ms 

Number of symbols K 

14 

1 

1 

Time-freq. product TF 

1.073 

1.25 

1.25 
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Fig. 12. Power spectral density. The power of the PRACH is chosen to achieve a similar PSD as PUSCH. 

B. Data transmission in PRACH 

Naturally, using the guard bands for data transmission causes an increased interference level in PUSCH. In Figure 
[T3] we show the effect on PUSCH symbol error rate caused by data transmission on a variable number of D-PRACH 
subcarriers, given the standard LTE PRACH and the new BFDM-based PRACH approach. 

Clearly, the performance of PUSCH does not deteriorate due to the proposed BFDM-based PRACH. By contrast, 
irrespective of the actual number of subcarriers used for data transmission, the BFDM-based approach leads to a 
slightly reduced interference level in PUSCH. Due to the strong influence of the D-PRACH on neighboring subcarriers 
in PUSCH, this effect is stronger when no DFT-spreading is used in PUSCH. The reason why larger gains, which could 
be expected from Figure cannot be realized is the PUSCH receiver procedure, which cuts out individual OFDM 
symbols from the received data. 
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Fig. 13. Symbol en'or rate in PUSCH (averaged over all (1200) subcaiTiers) plotted over the number of D-PRACH subcaniers. The BFDM-based 
approach slightly reduces the symbol error rate. This effect is stronger when no DFT spreading is performed in PUSCFI. 


C. Asynchronous users 

Asynchronous data transmission is a major challenge that comes with MTC and the loT. Therefore, we now consider 
a second, completely asynchronous, user that transmits data in the PRACH. Thereby we assign half of the subcarriers 
available for PRACH data transmission to this second user. However, we still evaluate only the performance of the 
original “user of interest” (and consequently we carry out channel estimation and decoding only for this user), which 
is assumed to transmit at the “inner” subcaiTiers close to the control PRACH. Thereby, we compare two waveforms, 
OFDM and the proposed spline approach. Figure shows the results. 

We observe that for completely asynchronous users, i.e., offsets larger than the CP (in which case OFDM loses its 
orthogonality property), the new pulse shaped approach reduced the symbol error rate up to a factor of almost one half. 
Nevertheless, the resulting symbol error rate may still seem excessive. However, as Figure [T^ shows, this effect can be 
compensated by allowing small guard bands (GB) in between the users. Figure [TS] compares the performance of no GB 
and GBs of up to 4 subcarriers, which already drastically reduces the symbol error rate. Interestingly, the spline-based 
approach without GB achieves roughly the same performance as OFDM with a GB of 4 subcarriers. In other words, 
we can save 4 subcatTiers using the spline-based PRACH. 

Although the symbol error rates depicted in Figure 14 and Figure may appear excessive, it should be noted that 


the delays considered in this evaluations (partially exceeding the cyclic prefix length) are unusually high. Figure 16 
gives a different picture, focusing on (positive and negative) delays around zero. 

































Fig. 14. Symbol error rate in PRACH (using 4-QAM) averaged over 10 out of 20 data subcarriers vs. the time offset of a second user. The other 
subcaniers are used by the second (asynchronous) user. The black line shows the CP length in LTE PRACH. 



Fig. 15. Symbol error rate in PRACH (using 4-QAM) over SNR with presence of a second asynchronous user. Here, the second user has a time 
offset of 200 /i.s. The spline based approach outperforms OFDM with or without a small number of guard bands. 
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Fig. 16. Symbol error rate in PRACH (using 4QAM) with perfect channel knowledge, averaged over 10 out of 20 data subcarriers vs. a varying 
time offset of a second user with a frequency offset of 62,5 Hz. The SNR is set to 25 dB. The black line shows the CP length in LTE PRACH. 


D. Simultaneous time and frequency offsets 

Figure [T6]illustrates the advantages of the spline-based transmission scheme for both frequency offsets and time delays. 
We plot the PRACH symbol error rates over the time offset of a second, asynchronous user. In addition, a constant 


small frequency offset of 62,5 Hz is applied. The SNR is fixed at 25 dB. Moreover, in Figure 16 we assume perfect 
channel knowledge for the user of interest. It can be observed that the additional frequency offset has a detrimental 
effect on both schemes, however, the performance loss of OFDM is significantly higher. 


VI. Conclusions 

We proposed and evaluated a novel pulse shaped random access scheme based on BFDM, which is especially suited in 
random access scenarios due to very long symbol lengths. It turns out, that the proposed approach is well suited to support 
data transmission within a 5G PRACH. In particular, numerical results indicate that the BFDM-based approach does not 
deteriorate PUSCH operations, in fact, it even leads to a slightly reduced interference in PUSCH when using (previously 
unused) guard bands for data transmission, irrespective of the number of subcarriers used for data transmission. Even 
more importantly, completely asynchronous users, with time offsets larger than the CP duration in standard PRACH, can 
be far better supported using the BFDM based approach than using standard OFDM/SCFDMA. The presented results 
will help to cope with the upcoming challenges of 5G wireless networks and the loT, such as sporadic traffic. 
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